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EXECUTIVE  SUMMARY 


The  goal  of  this  project  was  to  translate  measurements  of  perturbations  in 
parameters  of  sperm  motion  into  estimates  of  alterations  in  fertility  due  to  such 
changes  in  the  sperm.  These  sperm  changes  were  the  consequence  of  direct 
exposure  to  compounds  known  to  be  toxic  to  the  reproductive  function.  The  work 
was  organized  into  several  stages  and  was  carried  out  using  the  rabbit  as  a  model 
of  the  human  for  such  toxic  effects.  The  exposure  of  rabbit  spermatozoa  to  a 
series  of  compounds  was  performed  by  Dr.  Ronald  J.  Young,  U.S.  Army  Edgewood 
Research,  Development  and  Engineering  Center,  Toxicology  Branch.  Dr.  Young 
performed  replicate  experiments  in  which  videotapes  were  made  of  suspensions  of 
rabbit  sperm  exposed  to  a  range  of  concentrations  of  each  test  compound.  These 
tapes  were  sent  to  the  laboratory  at  the  University  of  California-Davis  (UC-D)  for 
subsequent  analysis  of  sperm  kinematic  parameters.  The  analysis  was  performed 
using  a  computer  vision  instrument,  CellTrak/S  (version  3.22,  Motion  Analysis 
Corporation,  Santa  Rosa,  CA);  the  generic  name  used  for  the  analysis  is  "CASA"- 
computer  aided  sperm  analysis.  An  initial  validation  study  was  performed  to 
evaluate  the  accuracy  of  CellTrak  measurements  for  rabbit  spermatozoa  as 
compared  to  manual  frame-by-frame  analysis.  After  the  system  was  validated  to 
UC-D  satisfaction,  each  compound  was  analyzed,  and  individual  files  were 
imported  into  Microsoft  Excel  (version  4.0}  for  data  manipulation  and  sorting.  The 
data  were  then  inputed  into  models  of  reproductive  risk  developed  at  UC-D,  again 
using  Microsoft  Excel. 
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PREFACE 


The  work  described  in  this  report  was  authorized  under  Contract  No. 
DAAA15-91-C-0130.  This  work  was  started  in  October  1992  and  completed  in 
March  1993. 

In  conducting  the  work  described  in  this  report,  the  investigators  adhered  to 
the  "Guide  for  the  Care  and  Use  of  Laboratory  Animals,”  National  Institute  of 
Health  Publication  No.  85-23,  1985,  as  promulgated  by  the  Committee  on  Revision 
of  the  Guide  for  Laboratory  Animal  Facilities  and  Care  of  the  Institute  of  Laboratory 
Animal  Resources,  Commission  of  Life  Sciences,  National  Research  Council 
(Washington,  DC).  These  investigations  were  also  performed  in  accordance  with 
the  requirements  of  AR  70-18,  "Laboratory  Animals,  Procurement,  Transportation, 
Use,  Care,  and  Public  Affairs,"  and  the  Laboratory  Animal  Use  and  Review 
Committee  (LAURC),  U.S.  Army  Chemical  Research,  Development  and  Engineering 
Center  (CRDEC),  *  which  oversees  the  use  of  laboratory  animals  by  reviewing  for 
approval  all  CRDEC  research  protocols  requiring  laboratory  animals. 

The  use  of  trade  names  or  manufacturers'  names  in  this  report  does  not 
constitute  an  official  endorsement  of  any  commercial  products.  This  report  may 
not  be  cited  for  purposes  of  advertisement. 

This  report  has  been  approved  for  release  to  the  public.  Registered  users 
should  request  additional  copies  from  the  Defense  Technical  Information  Center; 
unregistered  users  should  direct  such  requests  to  the  National  Technical 
Information  Service. 
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KINEMATIC  ANALYSIS  OF  RABBIT  SPERM  MOTION 
IN  THE  TOXICOLOGICAL  ASSESSMENT  OF  FERTILITY 

L  CASA  VALIDATION 

Purpose; 

To  evaluate  the  accuracy  of  CellTrak  measurements  for  VCL,  VSL  and  LIN  of  treated  and 
untreated  rabbit  sperm,  as  compared  to  a  manual  frame-by-frame  analysts. 

Methods; 

For  this  study,  the  Borax  Experiment  Tape  #2  (prepared  by  Dr.  Young  on  1/17/91  with  rabbit 
#339)  was  selected  from  the  available  tapes  due  to  the  optimal  quality  of  the  tape  with  regard  to 
optics,  sperm  concentration,  and  range  of  percent  motilities  for  the  treatment  dosages. 

A  tone-marker  was  audio-dubbed  onto  the  tape  at  the  beginning  of  each  field  after  the  image 
was  stabilized  and  focused.  The  tape  was  then  analyzed  by  CellTrak  at  60  Hz  using  one  digitized 
edge  and  the  default-sized  window  turned  on  to  minimize  data  overflows.  Analysis  was  triggered  by 
the  tone-markers,  and  the  following  Calibration  Set-up  parameters  were  used  (values  in  ( )  were 
used  for  analysis  of  "plucked"  30  Hz  data); 


Frame  rate  (frames/sec)  :  60  (30) 

Duration  of  data  capture  (frames)  :  30(15) 

Minimum  path  length  (frames)  30  (15) 

Minimum  motile  speed  (jj m/sec)  :  20 

Maximum  burst  speed  (pm/sec)  :  500 

Distance  scale  factor  (pm/pixel)  :  1 .83 

Camera  aspect  ratio  ;  1 .0000 

ALH  path  smoothing  factor  (frames)  1 1  (5) 

Cent.  X  search  neighborhood  (pixels)  :  4 

Cent.  Y  search  neighborhood  (pixels)  2 

Cent,  cell  size  minimum  (pixels)  1 

Cent,  cell  size  maximum  (pixels)  1 2 

Path  max.  interpolation  (frames)  :  1 

Path  prediction  percentage  (percent)  0 

Depth  of  sample  (pm)  :  10 

Video  processor  model  VP1 10 


A  hardcopy  was  made  of  the  sperm  tracks  in  each  field,  and  all  path  and  data  files  were  saved 
(test.pat,  vla.now,  mad.now).  Data  were  collected  for  25-50  sperm  for  each  of  the  following 
treatments:  Control  0.5  hr.,  Low  dose  0.5  hr.,  Medium  dose  0.5  hr.,  and  High  dose  0.5  hr. 

Manual  analysis  was  performed  on  the  same  sperm  at  approximately  the  same  time  as  were 
analyzed  by  CellTrak.  By  using  the  hardcopy  of  each  field  from  the  CellTrak  analysis  to  identify 
sperm,  and  the  audio  tone-marker  for  timing,  manual  sperm  tracking  began  within  several  frames  of 
the  beginning  of  the  CellTrak  analysis.  Using  a  16"  B/W  Electrohome  video  monitor  and  a 
Panasonic  AG6300,  the  head-midpiece  junction  of  each  sperm  was  tracked  frame-by-frame  and 
transferred  to  an  acetate  overlay.  Each  position  was  then  connected  by  a  line  to  create  a  path 
comprised  of  1 5  points  (30  Hz).  These  paths  were  then  measured  using  a  Numonics  electronic 
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digitizer  coupled  to  an  Apple  II  microcomputer  with  customized  software  that  calculates  VCL,  VSL 
and  LIN  for  each  sperm  track. 

The  scale-factor  for  the  manual  analysis  was  computed  by  measuring  the  image  of  a  Makler 
chamber  grid  on  the  monitor: 


Scale  factors  were  calculated  for  each  direction;  x,  y  and  diagonal: 

X  =  17.5  pm/cm 
Y*  15.7  pm/cm 
D  =  16.6  pm/cm 


Average  =  1 6.6  pm/cm 

Initially,  the  average  of  the  three  factors  was  used  for  the  manual  analysis.  (Previous  scale 
factor  calculated  for  Rat  Sperm  Validation  Study  was  16.7  pm/cm,  determined  in  the  X  direction 
only). 

Results: 

Examples  of  sperm  trajectories  for  all  treatments  and  both  framing  rates  are  shown  in  Figures 
1  (60  Hz)  and  2  (30  Hz).  The  control  and  the  low  dose  sperm  swam  with  very  linear  trajectories, 
while  the  medium  and  high  dose  sperm  swam  with  decreasing  linearities  and  increasing  ALH. 

Manual  vs.  CellTrak Analysis: 

Summary  statistics  and  correlations  for  each  parameter  and  all  treatment  groups  are  shown  in 
Table  1 .  Fifty  sperm  were  analyzed  for  the  control  and  the  low  dose,  while  only  33  and  25  sperm 
were  able  to  be  analyzed  at  the  medium  and  high  doses,  respectively,  due  to  the  lower  percent 
motilities  of  those  samples.  CellTrak  data  shown  are  for  the  30  Hz  analysis,  to  correspond  to  the 
manual  analyis  at  30  Hz.  Manual  values  were  calculated  with  two  different  scale  factors;  16.6  pm/cm 
(average  of  three  scale  factors,  MANUAL-1),  and  17.5  pm/cm  (x-direction  scale  factor,  MANUAL-2). 
Results  were  compared  between  manual  and  CellTrak  values  on  a  per-sperm  basis  using  a  paired-t 
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test;  all  parameters  were  significantly  different  between  Manual-1  and  CelfTrak  for  all  treatments  (p  < 
.01 ).  When  the  x-direction  scale  factor  was  used  to  compute  VCL  and  VSL  (MANUAL-2),  the 
medium  and  high  dose  VSL's  were  not  significantly  different  (p  >  .05),  and  the  control  VSL’s 
approached  non-significance  at  p  » .021.  LIN,  being  a  ratio  of  VSL/VCL,  remained  the  same 
regardless  of  the  scale  factor.  When  the  scale-factor  was  adusted  on  a  per-sperm  basis  according 
to  the  primary  direction  of  the  sperm's  trajectory,  all  manual  values  were  still  significantly  different 
from  CellTrak’s,  and  in  some  cases  were  even  more  significant. 


Manual-1 

CellTrak 

■uSSEXajl 

Manual-2 

CONTROL 

VCL 

108.914.6 

115.014.3 

.97 

114.814.9* 

.97 

VSL 

101.915.0 

103.414.9* 

.98 

107.415.2 

.98 

LIN 

.93 1 .02 

.88  1 .02 

.90 

(.93 1 .02) 

(.90) 

LOW 

VCL 

109.214.1 

116.014.0 

.96 

115.1 14.3* 

.96 

VSL 

101.914.4 

103.614.6* 

.99 

107.414.6 

.99 

LIN 

.92  1 .01 

.88  1 .02 

.68 

(.921.01) 

(.68) 

MEDIUM 

VCL 

71.216.4 

79.6 1 5.8 

.93 

75.1 1 6.8* 

.93 

VSL 

52.3  ±  5.6 

53.0  1 5.9* 

.97 

55.2 1 6.0* 

.97 

LIN 

.76  1 .04 

.62  1 .04 

.41 

(.76  1 .04) 

(.41) 

HIGH 

VCL 

77.2  1  5.8 

97.4  i  5.6 

.77 

81.416.1 

.77 

VSL 

51.315.1 

52.9 1 5.0* 

.96 

54.1  1  5.4* 

.96 

LIN 

.68  1 .05 

.56  1 .05 

.78 

(.68  1 .05) 

(78) 

Table  1.  Summary  statistics  and  correlation  coefficients  for  manual  and  CellTrak  values  of  VCL, 
VSL,  and  LIN.  Values  shown  are  mean  ±  SE.  Manual-1  values  were  calculated  with  the  average 
scale  factor,  16.6  pm/cm.  Manual-2  values  were  calculated  with  the  x-direction  scale  factor,  17.5 
pm/cm.  Numbers  in  bold  and  marked  by  *  =  nfll  significantly  different  at  p  >  .05. 

Regression  plots  for  manual  vs.  CellTrak  analyses  for  VCL,  VSL  and  LIN  are  shown  in  Figures 
3  (control  sperm),  4  (low  dose  sperm),  5  (medium  dose  sperm)  and  6  (high  dose  sperm).  The  high 
degrees  of  correlation  between  manual  and  CellTrak  datapoints  can  be  seen  in  these  plots  for  VCL 
and  VSL,  however  LIN  values  have  much  lower  correlation.  CellTrak  VCL  values  are  significantly 
higher  than  those  obtained  by  manual  analysis;  this  contributes  to  the  lower  LIN  values  reported  by 
CellTrak,  however  the  differences  between  manual  and  CellTrak  LIN  values  are  not  systematic. 
There  are  cases  where  the  manual  LIN  is  lower  than  the  CellTrak  LIN  for  the  same  sperm.  By 
careful  examination  of  the  individual  data  and  sperm  tracks  for  several  cases  where  there  was  a 
large  discrepancy  between  manual  and  CellTrak  LIN  values,  two  points  became  apparent.  First,  the 
machine  appears  to  be  much  better  able  to  discriminate  very  small  differences  in  the  sperm’s  head 
position  than  is  discernible  by  the  manual  method  of  tracking  the  head-midpiece  junction  (Figure  7). 
Second,  the  computer  is  better  able  to  track  erratically  moving  sperm  (Figure  8)  because  it  has  the 
advantage  of  detecting  the  head  position  in  relation  to  the  time,  whereas  by  the  manual  method 
employed  here,  the  tracking  (transferring  images  to  transparencies)  and  the  digitization  of  those 
trackpoints  are  two  separate  tasks,  usually  occuring  on  different  days.  Solutions  to  these  problems 
might  be  to  videotape  the  sperm  at  higher  magnifications,  so  that  minute  differences  in  head  position 
would  be  more  easily  seen,  and  also  to  utilize  a  digitizing  pen  to  track  the  sperm  directly  off  the 
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monitor  (Suarez  et  al,  1991 ),  rather  than  manually  drawing  the  tracks  and  measuring  them  at  a  later 
time. 

Effect  olframioa  Bale: 

Summary  statistics  for  each  parameter  at  60  Hz  and  30  Hz  are  shown  in  Table  2.  A  paired  t* 
test  was  used  to  compare  the  results  between  frame  rates.  All  parameters  except  VSL  were 
consistently  significantly  different;  VSL  was  significantly  different  between  the  two  frame  rates  for 
the  low  dose  sperm  only  (p  « .006).  This  difference  is  not  understood  and  cannot  be  explained. 
Regression  plots  of  the  data  (not  shown)  appear  to  be  similar  for  all  treatments,  yet  only  the  low 
dose  sperm  showed  a  significant  difference  between  60  and  30  Hz  values. 


FRAME  RATE 

VCL 

VSL 

LIN 

ALH 

CONTROL 

60 

141.9  ±3.6 

103.4  ±4.9 

.72  ±  .02 

5.78  ±.19 

30 

11 5.0  ±4.3 

103.4  ±4.9* 

.88  ±  .02 

4.18  ±.19 

LOW 

60 

142.9  ±4.3 

102.5  ±4.5 

.71  ±  .02 

5.77  ±.18 

30 

116.0  ±4.0 

103.6  ±4.6 

.88  ±  .02 

4.10  ±.19 

MEDIUM 

60 

129.8  ±6.3 

53.2  ±  5.9 

.39  ±  .04 

6.33  ±  .29 

30 

79.6  ±  5.8 

53.0  ±  5.9* 

.62  ±  .04 

4.32  ±  .33 

HIGH 

60 

141.3  ±6.0 

52.6  ±5.1 

.37  ±  .03 

7.08  ±  .45 

30 

97.4  ±  5.6 

52.9  ±  5.0* 

.56  ±  .05 

5.66  ±  .44 

Table  2.  Summary  statistics  for  CeliTrak  parameters  at  60  and  30  Hz.  Values  shown  are  mean  ± 
SE.  *  ■  not  significantly  different  at  p  >  .05. 

Table  3  shows  the  correlation  coefficients  between  the  60  and  30  Hz  data  sets. 


VCL 

VSL 

LIN 

ALH 

CONTROL 

.87 

1.0 

.91 

.57 

LOW 

.81 

1.0 

.92 

.44 

MEOIUM 

.90 

1.0 

.91 

.82 

HIGH 

.92 

.99 

.91 

.87 

Table  3.  Correlation  coefficients  for  60  vs.  30  Hz  CeliTrak  data. 


Conclusions: 

Rabbit  sperm  are  easily  digitizable  by  the  CeliTrak  instrument;  however  due  to  their  rapid, 
linear  motion,  care  must  be  taken  to  dilute  specimens  to  an  acceptable  concentration  range  in  order 
to  avoid  "losing"  sperm  to  collisions.  In  spite  of  the  fact  that  the  differences  between  manual  and 
CeliTrak  values  are  significant  in  most  cases,  especially  for  VCL  and  LIN,  there  is  generally  good 
agreement  between  the  two  methods,  given  the  limitation  of  the  manual  method  of  analysis  as 
described  earlier. 
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Text  continues  on  page  21 . 


FIGURE  1 

REPRESENTATIVE  SPERM  TRACKS  COLLECTED  AT  60  Hz 
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FIGURE  7 
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NOTE:  Manual  and  CellTrak  paths  are  not  at  the  same 
scale!!! 


H.  CASA  METHODS 


I 


CASA  was  performed  using  the  CellTrak  system  (Motion  Analysis  Corp.,  Santa  Rosa,  CA). 
Sperm  images  were  digitized  using  a  single  right  edge  (threshold  set  by  visually  matching  the  size  of 
the  digitized  image,  with  all  four  edges  on,  to  the  size  of  videotaped  images  of  progressively  motile 
sperm  heads).  Sperm  were  tracked  at  60  frames/sec  for  305  frames  (or  1/2  second).  Other  CellTrak 
calibration  setup  parameters  were  the  same  as  described  earlier  for  the  validation  study  without  the 
"window”  turned  on  to  allow  the  maximum  number  of  sperm  to  be  tracked  per  field. 

All  videotaped  fields  were  analyzed  for  each  treatment  in  order  to  maximize  the  number  of 
sperm  tracked.  However,  due  to  the  nature  of  the  experiments,  the  number  of  motile  sperm  varied 
according  to  dosage,  and  therefore  the  total  number  of  motile  sperm  analyzed  varied  as  well.  The 
number  of  motile  sperm  in  each  analysis  ranged  from  0  to  approximately  570,  and  averaged 
approximately  200.  Parameters  measured  were  %  motility,  straight-line  velocity  (VSL),  curvilinear 
velocity  (VCL),  linearity  (LIN),  amplitude  of  lateral  head  displacement  (ALH),  and  mean  angular 
deviation  (MAD). 

Following  is  a  list  of  each  compound  tested  and  videotaped.  Not  all  compounds  or  timepoints 
were  analyzed  by  CASA.  In  general,  only  compounds  where  three  experiments  were  repeated  at  the 
same  dosage  were  analyzed,  at  .5  and  1  hour  timepoints.  For  several  compounds,  a  2  hour 
timepoint  was  also  analyzed  if  no  effect  was  noted  at  the  earlier  timepoints. 


COMPOUND 

EXPT.# 

DATE 

RABBIT# 

DOSAGES 

CASA 

Propranolol 

1 

10/10/90 

465 

.5..1..01  mM 

X 

2 

12/07/90 

485 

.4,.  1  ,.05  mM 

(X) 

3 

02/05/91 

465 

.3,.  1..03  mM 

X* 

4 

04/28/92 

210 

.0003..0001, .00001M 

X* 

5 

04/30/92 

772 

.0003..0001..00001M 

X* 

Borax 

1 

12/10/90 

- 

32.5,13.0,6.5  mM 

not  sulvad 

2 

01/17/91 

339 

26.2,13.1,6.55  mM 

X* 

3 

11/18/91 

688 

26.2,13.1,6.55  mM 

X* 

Epichlorohydrin 

1 

01/29/91 

263 

12.7,6.35,3.175  mM 

X* 

2 

02/12/91 

465 

12.7,6.35,3.175  mM 

X* 

3 

02/13/91 

556 

12.7,6.35,3.175  mM 

X* 

PERC 

l 

01/14/91 

556 

3.91,1.96.0.978  mM 

notsulvad 
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<  .5  X 


Nonoxynol-9 


Chlorhexidine 


AlDha-Chlorohvdrin 


2,5-Hexanedione 


Dinitrobenzene 


EGME 


1 1/05/91 

712 

.00081, .0000405..0000162  M* 

11/07/91 

793 

.00081, .0000405..0000162  M 

11/12/91 

779 

.00081, .0000405..0000162  M 

04/16/92 

038 

.00081, .0000405,. 0000 162  M 

04/22/92 

779 

.00081, .0000405..0000162  M 

08/11/92 

772 

.00081, .0000405,. 00001 62  M 

10/17/91 

688 

.0236,.0ll8,.0059mM 

10/22/91 

793 

.0147,  .0059..00295  mM 

10/31/91 

854 

.00001 1 8,.0000059,.00000295  M 

05/04/92 

967 

.00001 18, .0000059, .00000295  M 

05/06/92 

972 

.00001 1 8, .0000059,. 00000295  M 

01/02/92 

793 

59.8,29.9,1 1.96  mM 

01/07/92 

712 

29.9,1 1.96,5.98  mM 

01/09/92 

772 

29.9,1 1.9,5.98  mM 

01/13/92 

688 

12.8,8.52,4.26  mM 

01/23/92 

038 

8.52,4.26,1.704  mM 

01/28/92 

779 

8.52,4.26,1.704  mM 

01/30/92 

712 

3.567,1.78,1.19  mM 

02/03/92 

213 

3.567,1.78.1.19  mM 

02/10/92 

793 

3.567,1.78.1.19  mM 

05/27/92 

772 

3.56.1.78,1.189  mM 

06/09/92 

779 

3.56.1.78,1.189  mM 

02/13/92 

779 

0.634,0.317.0.127  M 

02/18/92 

210 

0.634,0.317.0.127  M 

02/20/92 

712 

0.634,0.317,0.127  M 

06/1 1/92 

038 

.000025, .0000125..000005  M 

06/18/92 

712 

.000025..0000 125..000005  M 

07/16/92 

772 

.000025, .0000 1 25..000005  M 

07/21/92 

864 

.000025..00001 25..000005  M 

07/28/92 

793 

.000025..00001 25..000005  M 

X* 


X* 


X* 


X* 


X* 


X* 


notauKad 


not  niNad 


X* 


X 


X* 


notauhad 


X 


X 


not  aulvad 


X* 


X* 


X 


II  samples  were  taped  at  0.5, 1, 2, 4  and  6  hours  after  addition  of  the  compound.  Control  samples  were  diluted 
T6  medium  and  were  also  taped  at  t=0  hours. 

=  CASA  complete  for  .5  and  1  hour  (and  2  hour  for  propranolol  and  epichlorohydrin;  all  timepoints  for  lead); 
(X)  »  need  to  reanalyze  to  increase  sperm  numbers.  *  =  experiments  used  for  fertility  models. 
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After  each  analysis,  motile.ind and  motile.mad  ASCII  files  were  saved  and  stored  on  floppy  disks  for 
permanent  records  and  further  statistical  analysis.  These  files  contain  VCL,  LIN,  ALH,  and  MAD 
values  for  each  motile  (>20  pm/sec  VSL)  sperm  analyzed.  VSL  is  derived  by  multiplying  VCL  x  LIN. 
Test.pat  (binary)  was  also  saved  for  a  permanent  record  of  the  paths. 
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HI.  THE  RISK  ASSESSMENT  MODEL 


The  reproductive  risk  model  derives  from  ideas  initially  proposed  by  C.  van  Ouijn  for  predicting 
the  fertility  of  cattle  and  sheep  after  artificial  insemination  (van  Duijn,  1965).  van  Duijn  developed  a 
simple  model  relating  insemination  dose  to  the  probability  of  fertility.  We  have  extended  his  ideas  in  two 
basic  ways,  and  then  recast  them  into  a  model  of  risk  assessment  rather  than  fertility  itself.  First,  we 
have  utilized  a  body  of  knowledge  related  to  sperm  motility  and  transport  in  the  female,  and  to 
fertilization  itself,  that  was  obtained  after  van  Duijn's  pioneering  efforts  (cf.  Katz  et  al.  1989).  On  the 
basis  of  these  studies,  we  have  modified  the  criteria  for  functional  spermatozoa  to  include  thresholds  in 
measures  of  the  vigor  and  pattern  of  sperm  motion.  Thus,  our  “insemination  dose“  is  the  number  of 
sperm  that  have  a  velocity  >  100  pm/sec  and  a  linearity  of  motion  >  0.8  (see  below).  We  have  also 
utilized  more  modern  data  on  minimum  insemination  doses  necessary  to  achieve  fertility.  Here,  we 
considered  a  range  of  values,  in  order  to  account  for  the  difference  between  a  species  of  high  relative 
fertility  (the  rabbit)  and  a  species  with  lower  relative  fertility  (the  human).  We  have  retained  the  value  of 
the  parameter  kf  /  kq  =  1 .8,  obtained  experimentally  by  van  Duijn  for  cattle  and  sheep.  This  parameter  is 
the  ratio  of  the  temporal  decline  in  egg  fertilizability  divided  by  the  rate  of  decline  of  sperm  fertility  with 
time  (for  a  given  number  of  sperm).  We  believe,  on  the  basis  of  our  own  and  others'  studies,  that  this 
ratio  is  very  nearly  conserved  among  mammals.  We  did  conduct  parametric  studies  of  our  data  which 
demonstrated  a  relative  lack  of  sensitivity  to  small  perturbations  in  the  exact  value  of  this  parameter.  Our 
model  of  reproductive  risk  computes  the  ratio  of  the  altered  probability  of  fertility  (due  to  exposure  to  a 
toxicant)  divided  by  the  initial,  unperturbed  fertility.  This  ratio  is  akin  to  the  odds  ratio  frequently  used  in 
risk  assessment,  but  is  not  identical  to  it. 


1 )  For  a  fixed  delivery  flux  of  sperm  in  the  female  tract  to  the  egg: 


dF 

--^F  -F) 


or 


(1) 

(2) 
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F  •  probability  of  fertilization 

F‘  m  maximum  possible  value  of  F\F‘  <>  l 
t  -  time 

kf  m  rate  constant  for  change  in  fertilization  probability  with  time  for 
fixed  supply  of  sperm 

2)  Let  Q  =  Nv  be  the  net  transport  flux  or  supply  rate  of  fertile  sperm  in  the  female  tract 

N  s  number  of  fertile  sperm 
v  *  average  straight-line  velocity  of  fertile  sperm 


or 


(3) 


(4) 


Q,  a  initial  flux  of  functional  sperm 

ka  =  rate  constant  for  decline  of  sperm  flux  with  time 

3)  We  assume  that  there  is  a  limiting  or  threshold  value  of  Q,  call  it  Ox.  at  or  below  which 
there  is  no  chance  of  fertilization,  i.e. 

FaOand  QZQj  (5) 

Now  combine  equations  (1)  and  (3): 


dF  _kF  (F‘  -  F) 
dQ  ka  Q 

Integrating,  using  the  boundary  condition  (5) 


or 


F 


F' 


(6) 


(7) 
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We  shall  apply  equation  (7)  with  respect  to  the  initial  value  of  Q,  call  it  Q,,  at  insemination. 


Note  that  the  value  can  be  obtained  from  experimental  data  for  a  species  under  given 


conditions  of  insemination.  So  can  the  value 


°'X- 


Now  suppose  there  is  a  perturbation  in  sperm  numbers  and  quality  due  to  exposure  to  some 
toxicant.  Let  a  prime  symbol '  refer  to  these  conditions,  i.e 


Combining  equations  (7)  and  (8): 


F  1_{a 


i-<£)% 


is  the  reduction  in  fertility  due  to  exposure. 


and  kf/(  are  the  empirical  constraints  for  the  species. 

Q.  'o 


—  can  be  determined  experimentally  from  exposure  data. 


van  Duijn  calculated  the  values  -  1.7  - 1.8  for  different  data  sets  for  the  bovine.  These 
values  imply  that  the  decline  in  fertilization  probability  is  1.7  - 1 .8  X  the  decline  in  sperm  flux.  We 


still  have  to  input  the  ratio  for  the  species. 


Based  on  limited  data  in  the  human. 


1  Qr  1 
—  <  —  <  — 
20  Q,  10 
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For  species  with  higher  fertility  rates  (such  as  the  rabbit),  &  -  -~ 


For  our  exposure  data  for  the  rabbit,  we  compute  Q,  as  the  percentage  of  sperm  with  VSL  > 
100  pm/sec  and  UN  >  0.8,  multiplied  by  v,  the  average  value  of  VSL  for  this  subpopulation.  Thus,  for 
each  dosage  value  for  a  compound  we  compute  Qj.  Our  control  value  Q,  is  also  obtained.  These 

are  substituted  into  equation  (9)  to  determine  — . 

F 
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IV.  APPLICATION  OF  THE  RISK  ASSESSMENT  MODEL 

in  order  to  apply  the  fertility  model  described  above  to  our  experimental  data,  motHe.ind  files 
were  imported  into  Excel,  where  the  percent  of  progressive  sperm  (VSL  >  100  pm/sec  and  LIN  >  .8) 
was  calculated  for  each  sample.  The  model  was  applied  by  two  different  methods: 


Method^ 

For  all  compounds,  ~  was  calculated  for  each  treatment  (low,  medium  and  high  doses)  and 
experiments  as  follows: 


a  _(Pc)(VSLc) 

Q.  "  {Pr)(VSLr) 

Pc  =  %  progressive  for  Control  group 

VSLc  =  VSL  of  progressive  sperm  in  the  Control  group 

PT  »  %  progressive  for  Treated  group 

VSLf  =  VSL  for  progressive  sperm  in  the  Treated  group 


Then,  using  ratios  of  —  (or  .05)  and  —  (or  .01),  and*^/  =  1.8,  —  was  calculated  for 
Q.  20  100  /*b  F 


Q. 

each  treatment  from  formula  (9): 


F_ 

F 


1- 


(g-)%](§-)% 


l-(^) 


% 


Results  of  —  calculated  by  this  method  for  each  compound  follow. 
F 
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Text  continues  on  page  42. 
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ALPHACHLOROHYDRIN  DOSE  RESPONSE  OF  F/F  FOR  Q»/Qo  *  1/20  AND  1/100 
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DOSE  RESPONSE  OF  F/F  FOR  Qt/Qo  =  1/20  AND  1/100 


CHLORHEXID1NE  DOSE  RESPONSE  OF  F  /F  FOR  Qt/Qo  ■  1  /20  AND  1/100 


DNBFERT.XLS 

DINITROBENZENE  DOSE  RESPONSE  OF  F/F  FOR  Qt/Qo  =  1/20  AND  1/100 


DOSE  RESPONSE  OF  F/F  FOR  Qt/Qo  =  1/20  AND  1/100 


EPIFERT.XLS 

EPICHLOROHYDRIN  DOSE  RESPONSE  OF  F/F  FOR  Qf/Qo  =  1/20  AND  1/100 
.5  HR 


If  38  38  88 

5«-  o  -  o  -  ~  - 


5*  88  88  8:8 

2fc  "  " 

HI 


^  §  8 5  85  8 5 

5  5  do  do  do 
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HEXAFERT.XLS 

2.5  -  HEXANEDIONE  DOSE  RESPONSE  OF  F/F  FOR  Qt/Qo  «  1/20  AND  1/100 


LEADFERT.XLS 

DOSE  RESPONSE  OF  F/F  FOR  Qt/Qo  =  1/20  AND  1/100 


88 

001 

001 

88 

88 

88 

§8 

88 

88 

88 

100 

soo 

85 
o  6 

8  o 
o  o 

*2 

2 

-0  o 
S  £ 

45% 

121.6 

59% 

138.4 

ti 

-  cO 

g  jo 

a#  ^ 

ex  2 

SI 

K  « 

.  «? 

28% 

1264 

£§ 

:  88 

88 

88 

88 

88 

88 

§8 

88 

88 

100 

soo 

100 

soo 

0.05 

0.01 

51% 

135.9 

Si 

20% 

112.3 

35% 

135 

44% 

132.9 

30% 

123.6 

40% 

130 

44% 

131.9 

27% 

122.6 

20% 

125.6 

38% 

125.5 

22% 

1126 

S  n 
2  »|  * 


—i  s 


36 


IEADFERT.XLS 
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NONOXYNOL-9  nofertxls 


PROPFERT.XLS 

PROPRANOLOL  DOSE  RESPONSE  OF  F*/F  FOR  Qt/Qo  =  1/20  AND  1/100 


122.3  I  1157  I  118.5  1  1 10.3  |  0.01  |  1,QQ 


Method-2 


In  order  to  eliminate  some  of  the  variation  between  experiments  seen  when  ~  was  calculated 
by  Method  1 .  a  second  method  was  employed.  This  time  was  calculated  for  each  treatment  as 

F 

described  for  Method  1  and  then  averaged  across  all  experiments.  Then,  a  single  value  of  —  was 

F 


calculated  from  formula  (9}  lor  each  treatment,  using  &-s-L  only.  Results  using  this  method  follow. 

Q,  20 
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Text  continues  on  page  68. 
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DOSE  RESPONSE  OF  F*/F  FOR  Qt/Qo  =  1/20 


DOSE  RESPONSE  OF  F'/F  FOR  Qt/Qo  =  1  /20 


DOSE  RESPONSE  OF  P/F  FOR  Qt/Qo  =  1/20 


DOSE  RESPONSE  OF  F/F  FOR  Qt/Qo  =  1  /20 
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EGME1MDL.XLS 


EPICHLOROHYDRIN 
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EPICHLOROHYDRIN 
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DOSE  RESPONSE  OF  F'/F  FOR  Qt/Qo  =  1/20 
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DOSE  RESPONSE  OF  P/F  FOR  Qt/Qo  =  1/20 
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DOSE  RESPONSE  OF  F*/F  FOR  Qt/Qo  =  1/20 
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PROPRANOLOL 

NOTE:  DOSAGES  DIFFER  FOR  EACH  EXPERIMENT!! 

EXPT.  1  »  EXPT.  3  -  .3..1..03mM;  EXPT.  4.5  -  .3,.1..01mM 


DOSE  RESPONSE  OF  P/F  FOR  Qt/Qo  =  1/20 


DOSE  RESPONSE  OF  F  /F  FOR  Qt/Qo  =  1/20 
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V.  C0NCLU8I0N8 

The  compounds  studied  exhibited  a  range  of  effects  upon  sperm  motion  and  subsequent  fertility. 
Our  selection  of  dose  levels  was  designed  to  span  a  range  which  just  attained  a  threshold  In  visually 
obvious  diminution  of  sperm  motility.  That  is,  we  sought  to  study  the  onset  of  deleterious  effects  on 
sperm  and  fertility.  The  objective  analysis  of  sperm  motion  and  risk  assessment  produced  a  range  of 
effects  for  the  different  toxic  compounds.  This  range  exemplifies  the  importance  of  such  objective 
analysis,  as  compared  to  visual  estimates  of  cessation  of  sperm  motility;  according  to  the  latter,  the 
doses  for  all  compounds  produced  relatively  similar  effects.  Given  that  the  concentration  ranges  of  the 
different  compounds  applied  to  sperm  varied,  we  can  rank  their  level  of  potency  as  follows: 

Mighty  potent  Dinitrobenzene 

Nonoxynol*9 

Potent  Borax 

Propranolol 

EQME 

Moderately  potent  Alphachlorohydrin 

Eplchlorohydrin 

2.5-Hexanedione 

Lead 

This  ranking  takes  into  account  differences  in  the  time  course  of  action  against  the  spermatozoa,  as  we 
as  the  dose  responses  for  each  time  point.  For  example,  lead  exhibited  virtually  no  effect  until  extender 
Incubation  with  sperm. 

It  should  be  appreciated  that  these  results  are  for  an  In  vitro  bioassay  of  effects  of  compounds 
applied  topically  to  spermatozoa.  Such  effects  may  not  be  identical  to  those  that  occur  when  the  entire 
organism  is  exposed  to  such  compounds.  In  addition,  effects  on  the  human  (in  vitro  as  well  as  in  vivo) 
may  not  be  the  same  as  those  in  the  rabbit.  However,  the  rabbit  is  a  commonly  used  model  in  testing 
compounds  that  are  chemical  contraceptives.  Thus,  its  use  in  this  context  is  reasonable.  Overall,  this 
study  has  presented  a  new,  objective  approach  for  using  sperm  as  biomarkers  of  reproductive  risk 
assessment.  It  has  provided  estimates  of  the  onset  of  risk,  with  increasing  exposure,  that  would  not 
have  been  possible  with  traditional  visual  methods  of  sperm  assessment  and  statistical  analysis.  We 
hope  that  this  approach  will  be  of  value  in  further  analysis  of  reproductive  risk  after  toxic  exposure  of  a 
male,  e.g.  in  screening  potentially  deleterious  workplace  or  environmental  exposures. 


